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Abstract 

Decay amplitudes for mesons in two-dimensional QCD are discussed. We 
show that in spite of an infinite number of conserved charges, particle pro- 
duction is not entirely suppressed. This phenomenon is explained in terms 
of quantum corrections to the combined algebra of higher-conserved and 
spectrum-generating currents. We predict the qualitative form of particle 
production probabilities and verify that they are in agreement with nu- 
merical data. We also discuss four-dimensional self-dual Yang- Mills theory 
in the light of our results. 
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1 Introduction 



Major advances in two-dimensional QCD have been made since the 
pioneering work of 't Hooft M, who found the mesonic spectrum of the 
theory in the limit of a large number of colours. In subsequent works, an 
exact expression for the decay amplitudes of the mesons was found [^, ^] . 
However, both the 't Hooft equation and the decay amplitudes could only 
be solved numerically. The numerical solution of the theory has been 
enterprised by many authors j|, ||. The numerical calculations are difficult 
to converge for massless fermions. A strongly convergent algorithm for 
very reliable numerical computations is also available in the literature! 14]. 
These new results open challenging ways of testing various aspects of the 
theory, specifically those concerning decay amplitudes of the theory. 

It has been suggested [|J that massless two-dimensional QCD is an 
integrable system. It is believed that integrability of a theory implies 
stability of its bound states. One could, therefore, expect vanishing decay 
amplitudes for the mesonic states of the theory. The numerical results 
indicate otherwise. The decay amplitudes, however, are small relative to 
those obtained for the massive theory. 

The main purpose of this article is to resolve this apparent conflict by 
relating diverse properties of two-dimensional QCD, as well as understand 
more about the dynamical structure of two dimensional QCD. Specifically, 
we show that the 't Hooft sector and the integrable sector of the theory 
are decoupled, in the sense that the conservation laws found before do not 
commute with the spectrum generating algebra due to quantum correc- 
tions. We also believe that the theory is not integrable in the 't Hooft 
sector, that is to say that the integrability properties are spoiled by quan- 
tum corrections. However, we think that a strong simplification occurs 
in the massless case as compared to the massive fermion theory, allowing 
us to draw special attention to such a difference in terms of the conserva- 
tion laws and the anomalies, which presumably have a mild contribution 
semi-classically. 

We review the known bosonisation procedure in Section 2 and show 
how the theory can be re-expressed in a fermionic language, by using 
gauge-invariant chiral fermions. The final re-fermionised action provides a 
convenient starting point for the analysis of the mesonic spectrum. The 
significant outcome of this procedure is that the Poisson brackets receive 
quantum corrections. The conserved Noether currents then satisfy the 
Kac-Moody algebra only on taking the contribution from the anomalous 
Poisson bracket into account. In Section 3, we combine the algebra of 
these conserved currents with that of the fermionic bilinears, which have 
been shown to generate a Woo algebra. We show that, except for the pion, 
the mesons generated by the above W m algebra are not eigenstates of the 
Sugawara operator. Conservation of the Sugawara operator indicates the 
existence of an infinite number of conservation laws. If such conservation 
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laws survive quantization and physical states are eigenstates of the Sug- 
awara operator then the dynamics of the theory is severely constrained. On 
the other hand, the breakdown of conservation laws by anomalous terms 
permits the decay of higher states. This explains the small decay rate nu- 
merically observed in the massless theory. Based on these results, we make 
detailed predictions, in Section 4, about various features of the spectrum. 
In Section 5, we use the recent numerical results to verify our predictions 
and hence the accuracy of our calculational results. 



2 Bosonisation vs Fermionisation 



L = --F vaf F* a ' + virfD v ri>, (1) 



We start with the lagrangian 

1 

T 

which describes QCD with massless fermions. In two dimensions, the the- 
ory can be bosonised by computing the fermionic determinant which arises 
in performing the path integration. In the bosonisation procedure one 
first defines left and right components of the gauge field in terms of the 
matrix- valued fields U and V, i.e., 

A+ = -U- l d + U and A. = -Vd-V' 1 . (2) 

e e 

It is well-known that the Jacobian resulting from the above change of 
variables together with the effective action can be expressed in terms of 
the WZW functional @. That is to say, 

deti $ = j w[A] = Jv 9 e iSF ^. (3) 

The bosonic action Sf[A, g] may be obtained by the repeated application 
of the Polyakov-Wiegmann identity || [| and is given by, 

S F [A,g] = T[g} + ^-Jd 2 x[e 2 A 2 -e 2 A + gA-g- 1 

-ieA+gd-g- 1 - ieA-g^d+g] 
= T[UgV]-T[UV], (4) 

where T is the WZW action functional. Subsequent inclusion of the Fadeev- 
Popov ghost action S[gh], replacement of the variables UV by S and UgV 
by fusing systematically the invariance of the Haar measure, and finally 
the introduction of a scalar field E, to disentangle the F^F^ interaction, 
yields the following form of the partition function || 0] : 



V0 T ^] V ^ gh yS[gh} VTiVEe -i(c v +i)T[n 
exp [-(cy + l)tr / d^d+^Ed-E" 1 - 2e 2 (c v + l) 2 tr / d 2 xE 2 }. 



(5) 
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From a further replacement of the variable E by a field /?, satisfying the 
relation d+E = j^f3~ 1 d+/3, the factorised form of the partition function, 

Z = Z[g, l]Z[% -(c y + l)]Z gh Z p , (6) 

is obtained. In the above expression, Z[M,n] is the partition function of 
a WZW field M with central charge n, Z g ^ is the ghost contribution, the 
non-trivial coupling-constant-dependent part of the partition function is 



Z. 
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and the auxiliary field C_ is introduced to make the action local. []. The 
equations of motion are equivalent to the conservation law <9+i_ = 0, where 

= 47re 2 /r 1 d + /3 - 2U'_ + z(4vr) 2 e 3 C7_ - (4vre) 2 [C_, IL_] , (8) 

14 = <9 + C_ and Ilg = ^<9_|_/? _1 -\-ieC-f3~ 1 are canonical momenta conju- 
gate to C- and /? respectively. Canonical commutation relations lead to a 
Kac- Moody algebra obeyed by the current Z_. As in the case of conformally 
invariant theories, the conservation law <9+/_ = implies that there is an 
infinite number of conserved charges, = J I^(x~)(x~) n dx~ . These 

act on asymptotic fields as multiplications by powers of the (negative) 
component of the momentum times right-SU(n) transformation, i.e., 

Q (n) ij |0M(p)) nu p™5 jl U jn \(3 kn ). (9) 

This puts a stringent requirement on particle scattering, forbidding particle 
production or decay f\ 



Recently, there have been suggestions [14] that the stability of the spec- 
trum can be tested in the framework of 't Hooft's mesons decay amplitudes. 
The numerical computations of the amplitudes provide a detailed check of 
the claims presented here, namely that the infinite number of conservation 
laws does not imply absence of particle production or stability of the meson 
and zero decay amplitude. 

As mesons are fermionic bound states, the relevance of the above con- 
served quantities to the decay of mesons can be conveniently studied by 
writing the (5 sector of the theory, defined in ([?]), in terms of fermions. In 
the fermionisation procedure, the WZW term is equivalent to an action for 
massless fermions and the interactions are treated in perturbation theory 
by using the adiabatic principle of form invariance |H 



1 The factorised form of the partition function is actually elusive, due to several BRST 
constraints remaining from the gauge condition and change of variables. For the definition of 
vacuum and physical states the reader is referred to jl0[ [ll| [l2| . 

2 A detailed discussion of the action of conserved currents on (3 field, and the consequences 
for /^-scattering have been discussed in ref. |0, where an exact S-matrix has been conjectured. 
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Using these techniques, we obtain the following fermionised form of the 
action: 

S = J d 2 x[^i(5 ji d. - ie&i)^ + itr(3+ C_) 2 + i^ld+^i}. (10) 

The 't Hooft spectrum can be algebraically generated by using the 
spectrum-generating algebra of gauge-invariant fermion bilinears. Since 
the (3 fields in (0) , as well as their fermionic replacements in ([k]) , are gauge- 
invariant objects, the above action represents a chiral theory which has to 
be quantised using anomalous Poisson brackets |l6| (see also chapters 13 



and 14 of [15Q. Indeed, the left-moving Noether currents, 

it = i>H - -n'_ - -^c_ + i[c_,n_], (11) 

e 47r 

generate a Kac-Moody algebra only if the anomalous Poisson Bracket 
(APB), 

{n-(x),VV(y)} APB = ^(ar-y) J (12) 

is used J]. 

The operator (|ll]) (or (|8|)) contains colour indices. However, in 't 
Hooft's formulation one deals exclusively with colourless states and, there- 
fore, it is natural to consider the Sugawara operators 

L{x) =: IyWljiix) : , (13) 

where the colour indices are contracted and the divergent terms in the Wil- 
son expansion (corresponding to the Kac-Moody algebra) are subtracted 
by means of a normal-ordering prescription. 



3 Classification of the mesonic states 

In this section, we review the action of the algebra on the spec- 
trum. The spectrum-generating algebra, as obtained in |l7j, arises from 

the bilinears M a ^(x,y) = tpj Jl (x)e ie ^ Afl ^ d ^ip^(y), where a and 7 are the 
chirality indices (+,—). In the massive case, these bilinears are related by 
the equations of motion. In the massless case, however, the mixed term, 

M-| , corresponds to the integrable-/3 sector, while the right-right term, 

M , corresponds to the usual meson bound states of 't Hooft. Therefore, 

we take 

M__(:r~,y~;2; + ) = if>*l(x~, x + )^i(y"; x + ) (14) 

as the spectrum-generating current. In the momentum space of the x~ 
variables, a classical solution of the equations of motion obeying the Gauss 
constraint is 

(M__(fc_X ;x+)) class = 5(k^ - (15) 

3 The bosonic formulation is advantageous since it contains information of order Ti already 
at the classical level. 
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Expanding M 



around such a solution in the large N limit 

M___ = eVW w (M_.) class e-VW w 



17] and using 
(16) 



one finds that the Fourier modes of W, which represent quantum fluctua- 
tion around classical solution, obey the 't Hooft equation. Thus, the W m 
algebra found for the bilinears in [l?]] is a spectrum-generating algebra. 
We emphasise that in the massive theory, the full content of the theory 

is preserved in the individual chiral sectors, M++ and M , which are 

related to each other by the fermion equation of motion. In the massless 
case, however, we study the bilinears constructed from r ip + . There is a 
mixing of the sectors and it does not suffice to study one sector on its own. 



4 The combined algebra 

In this section, we study the interplay between the Woo-spectrum- 



generating algebra and the conserved current (13). We identify the bi- 



linear ip\_ip+ as the pion and ip^_D n i(j + as the higher states obtained from 
the fermion bilinear M ++ . 

The pion is an eigenstate of L(x) Q 

2 2 

[L{x),i>\y)i>(y)} = 6 ^ip(y)+ 6 M ^{y)) 

2{x - y) A 2{x - y) 

+ regular terms. (17) 

The higher-state fermion bilinears, however, due to various anomalous 
terms, are not eigenstates; 

[L(x),dip t ip(y) - V> f (y)dip(y)} = 
e 2 2 

;(dijjty(y) -ipWip(y)) 



4 (x — y) 2 



— — d(d^ip(y) — ip^dip(y)) 
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4 J 



x-y (x - y) 

+ e 2 L 1 , 9 (V>*WV) — (V*VWV))] 

(x-yy x-y 

+ [ - — [U kj C ki - C jk U ik ] - 16—C ji + edilF] x 

2 7T 

x — y [x — y) A 

-l-anomalous Poisson Brackets terms + regular terms. (18) 

Using the above two equations, we can write down the general form 
of the algebra obtained from the action of the Sugawara operator on the 
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Recall that in the present fcrmionic formulation we have to use the anomalous Poisson 



bracket. 
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physical states, that is, 



n + 1 1 

L ■ M n rs r^M n H 9M n + anomalous terms. (19) 

(x-yy x-y 

Thus we see that, in the absence of anomalies, the n" 1 state is an eigenstate 
of L with eigenvalue n+1. Taking this into account, we find that the action 
of L on the in and out states is given by 



(n - n\ - n 2 - l)(M m M n2 \M nj ) = J dx(\L out - L in \) 

/roc (I 
dx / dt — (M ni M n2 1 L | M n ) = anomalies, (20) 
J —oo v& 

which shows the existence of an infinite number of conservation equations. 
In the absence of the anomalous terms, the decay amplitudes vanish. How- 
ever, when anomalous terms are present unphysical (non-mesonic) oper- 
ators are introduced in the right-hand side. Therefore, the amplitude is 
non-vanishing. If one of the out states is a pion, then we will have the 
following set of recursive relations: 



(n-m-l){M M ni \M n ) = J2( M o X \ M n) 

x 

+ J2(M M ni \Y) (21) 

Y 

where X, Y are bound states of lower states and a pion state is always 
present on the right-hand side. The solution to the above relations is 
found by using 

(M M \Y) = 0, (22) 

which is obtained from the 1/N expansion of the decay amplitude and 
is valid up to second order. These solutions require the vanishing of the 
decays involving pions 0. 

From eqn. (|l8|), we see that there are further corrections which come 
from both higher terms in the Wilson expansion and anomalous Poisson 
bracket quantization of the chiral fermions. These terms cannot cancel 
one another, due to their different functional form (as an example, higher 
II derivatives never arise from anomalous Poisson Brackets). There are 
further higher terms when one considers the gauge-covariantized current. 

The decay amplitudes in the massless theory (which is integrable in the 
absence of anomalies) are suppressed as compared to those of the massive 
theory (non- integrable). This is entirely due to the quantum corrections. 
Since the anomalous terms are of order h, they disappear for quasi-classical 
decays. In the massive theory, in addition to quantum corrections, there 
are mass terms which spoil integrability. These terms do not vanish in the 



It is worth mentioning that the algebra is independent of the number of colours, and the 
anomalous term only makes a first-order contribution in the 1/N expansion. 
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quasi-classical approximation. In the massless theory, this approximation 
is reliable for decays of highly-placed states (e.g. near the continuum limit ) 
to large-momentum states. On the contrary, decays into small-momentum 
states, marked by interference terms, are highly influenced by quantum 
corrections which spoil quasi-classical approximation. 

On the basis of the preceding results, we predict the following features: 

1. The pion decouples in the massless fermion theory. Such a decou- 
pling is valid in all orders of the 1/N expansion. 

2. For massless fermions, the decay of large-mass states into states with 
large momenta is severely suppressed. 

3. The probabilities for further decays, although reduced, are not very 
small. For small N, the amplitudes are much smaller than the corre- 
sponding ones in the massive case, due to the enhanced significance of the 
larger-order terms in the 1/N expansion []. 

4. Finally, the decay amplitudes of very massive mesons vanish, unex- 
ceptionally, in the massless fermion case. 

5 Numerical back-up and large N behaviour 

In this section, we verify these predictions by means of the numerical 
computation of the amplitudes (for predictions 1,2 and 3), as well as, by 
using large ./V analysis (for prediction 4). 

Amplitudes for meson decay were initially derived in the framework of 
1/iV-expansion in [g], and in more detail, including higher order correc- 
tions, in [||. The 1/N corrections vanish in the massless case Q This can 
be explained by studying the expression for the decay amplitude ||, 

A = (i-c)— !- r axMx)M-)%(^) 

1 — UJ Jo UJ 1 — UJ 

1 X UJ X 

-(l-C)-dx dx<f>n(x)<l> q £ )$„(-) 

UJ Ju) 1 — UJ U 

1 f ( u T 

+-(1-Q-I2- / dxd> n (x)M-), (23) 
N 1 — uj Jo uj 

where uj = <fr n (x) = J dy ^yji and C denotes the interchange of final 
states. 

We observe that the higher-order corrections are always multiplied by 
the factor 

f n = [ dx<j) n {x), (24) 
Jo 

6 We recall again that 1 /N corrections are very important in the massive case which strength- 
ens these predictions for groups such as SU(2). 

7 The lowest-order term, however, seems to give a non-vanishing contribution. This is dis- 
cussed later in this section. 
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where <f> n {x) is 't Hooft's wave function of the decaying state. It can 
be verified using 't Hooft's wave equation, that f n vanishes for massless 
fermions 0]. Moreover, the authors of ref. |J claim that, as higher-order 
corrections amount to a redefinition of constants in the massless case, only 
the very first term in the 1/N series survives. Because this does not hold in 
the massive case, corrections to all order exist. We restrict our analysis to 
the large- N limit where amplitude for both massless and massive fermions 
may be non- vanishing. 

In order to fully understand the behaviour of the leading-order term 
(i.e. the first term in expression (p3|)) we need to solve the 't Hooft equation 
numerically. The solution is then inserted back into (^) and the numerical 
integration is carried out. We are then in a position to compare the results 
for the massless and the massive cases. 

The method which employs a MATHEMATICA program is unsatisfac- 
tory because of its bad-convergence behaviour. The numerical results are 
unreliable although they are compatible with the predicted vanishing decay 
rate in the massless case, and the non- vanishing decay rate in the massive 
case. Therefore, a more sophisticated method of numerical computation is 
needed ^ . 



A more elaborate method [14| confirms that particles do not decay 
into states which contain the pion whereas in the massive case, the decay 
amplitudes are shown to be significantly large (see Figure 1). 





Figure 1: Decay rate of the 5 th excited state into the 1 st and ground (pion) states 
is plotted as a function of the outgoing momenta. The first diagram corresponds 
to massless and the second to massive fermions. The vertical bars mark the on- 
mass-shell values. 



For the massless case, several decays of the type n — > 1 + 1, where 1 



8 For a survey of these methods using orthogonal eigenfunctions expansion we refer the reader 

to apira. 
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denotes the first-excited level of 't Hooft's series, have been computed (see 
Figure 3). For higher values of n the amplitude approaches zero rapidly. 
For the massive case, the amplitudes vary randomly (see Figure 4). 

The overall results can be summarised as follows : 

(i). In the massless case, for a general decay k — > I + p the amplitude 
can start from a high value (see specifically the example 14 — ► 3 + 3 in 
reference |]i~4j| ) but decreases rapidly with increasing k . These results 
are also compatible with the earlier results of reference ||] . Although the 
results for the massless case are very precise, those for the massive case 
require further refinement. However, all the results obtained are precise 
enough for the conclusions drawn hitherto, and will be verified by other 
methods. 

(%%). An exact computation of the decay rate of the higher-mass states 
guarantees the correctness of the numerical results obtained so far. Indeed, 
using the large N techniques of Q or [^0| one finds the required amplitude 
in the large k limit. In terms of the fermion mass m, a parameter 7 
satisfying 

irm 2 

7T7COt7T7 = 1 ^— - (25) 

is defined (the massless case corresponds to 7 = 0). In the large n limit, 
the decay amplitude is found to obey 

A ~ sin7T7, (26) 

which is valid up to constants describing the behaviour of the 't Hooft 



wave functions at the origin (see |2C] for further details) and which shows 



a non- vanishing result only for non-zero values of the parameter 7. The 
same computation may be explicitly carried out in the massless case. 

In order to investigate the issue of pair production for large number of 
colours and high mesonic states, we calculate the amplitude of the decay 
of a very massive meson of mass m n into two mesons of masses m ni and 
m n2 . Such an amplitude, at lowest order in the inverse number of colours 
and up to an overall constant, is given by the vertex 

V ni n 2 n = e 2 — f dx [ dy 1 

P2 JO JO (X^) 

x i> n M4>nM\^ L ]r v , )-u \ + i xp i P2 )}, (27) 
1+P1/P2 1+P1/T2 

where pi is the plus-component of the momentum. For decay of very 
massive mesons, m n — > 00, energy-momentum conservation implies 

ml = n^. (28) 

P2 

Thus the main contribution to the integral comes from x ~ 0. Next, we 
make the following substitutions 

* = (29) 
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(t>i(x) » C iX \ (30) 
where 7 is defined in (^) and x is small, and 

Mt(— ) 2 )«^)=sm(- + 5(6), (31) 
where <5(£) is a phase f]un equation (p7|) to find 

Vn^n = ( — ) 27 d r^/^xV 7 ^^[0((l-x')4)-^(( 1 + ^) ! 4)] 

= C^y^dh. (32) 

e 

Since the vertex is symmetric under the exchange 1 <-> 2, the relation 
Cj/Ij = r must be the same for i = 1, 2 and can be computed in a conve- 
nient limit. We thus compute r for i = 2 and large values of mi- In this 
limit, the asymptotic behaviour of the wave functions can be used to find 
the explicit values of C2 and I2, i.e., 

C 2 « ( ) , I2 « ( ) 'vrsmv^. (33) 

e e 

Subsequently, the expression 

V nin2 n ~ e 2 (— ) 27 sin tt7 (34) 

e 

is obtained which shows that the vertex vanishes for 7 — > 0. 

For massless fermions, that is for 7 = 0, the rj integration in the first 
term of eqn. ( |32"|) can be performed. Moreover, by using 't Hooft's equation 
in the same term (to replace | ) we obtain 

1 dxfaix) /°°dry ^, (1 + 7? ^ 2) - ( l dx4> 2 {x)[ii 2 + -^—](t>{{l-x)ii 2 ) 



(77 + x) 2 Jo 1 - x 



/* j0 (a; - 1 + -t) 



where /j = ^2.. The last two terms cancel [^lj due to the identity 

^7/% = -(- + 1)^)- (36) 

In order to confirm these conclusions further, we make a detailed com- 
parison between the amplitudes obtained in the massive and massless cases. 
This can be conveniently done by considering the table presented below. 



The validity of equation (j3lj) is confirmed for massless fermions by using numerical 
simulations. 
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Decay series 


A 


fc 


.A/Jfc 


8 -> 1 + 1, m = 


.25 


.4 


.65 


8 -» 1 + 1, m + 


.5 


.25 


2 


9 -> 1 + 1, m = 


.2 


.4 


.5 


9 -»■ 1 + 1, m + 


.5 


.28 


1.8 


10 -> 1 + 1, 


.8 


.3 


2.7 


11 -> 1 + 1, m = 


(.05 


.45 


(•1 


13 -»■ 1 + 1, m = 


sa 


.45 


« 



In this table, the second column represents the momenta of the outgoing 
particles and the figures in the last column are proportional to the decay 
probabilities 0. These tabulated results once again confirm the prediction 
that the ratio between massless and massive fermion decay rates goes to 
zero-being smaller than 3% for the 11th state. For other series of decays, 
this ratio approaches zero more slowly. Nevertheless, one clearly observes 
that this ratio approaches zero, e.g., for massless series k — > 2 + 2 (see F 
igure 2). 











V 


v v 




Figure 2: Decay rates of the 10 th , 11 th and 12 th states into two mesons in the 
2 nd excited state (for massless fermions) are plotted as functions of the outgoing 
momenta. The vertical bars mark the on-mass-shell values. 



10 We have chosen the units such that the mass of the initial state is unity (the influence of 
the fermion mass is small in the present cases). 
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6 Conclusion 



Massless QCD contains higher-conservation laws which in general imply 
integrability. These conservation laws have been derived in the massive 
(/?) sector. The mesons in the 't Hooft sector are built up of fermion 
bilinears which are dressed with bosonic fields of the massless sector. We 
have shown here that the spectrum-generating algebra, which defines the 
't Hooft sector, does not commute with the higher conservation laws due 
to quantum corrections to the short-distance expansions. This implies the 
break down of integrability in the meson sector. The quantum nature 
of these corrections means that they are insignificant for quasi-classical 
decays. This renders the theory quasi-integrable and accounts for the exact 
decoupling of the pion. 

The theory has a complex structure of constraints. The mesonic states, 
although physical, are not eigenstates of the conserved charges obtained 
from the Sugawara operators. This is because the massive and the mass- 
less sectors of the theory are connected by the constraint equations. The 
massive sector alone (the (3 sector) is integrable but does not generate 
the physical Hilbert space since it is not BRST invariant. The change of 
variables made to decouple the dynamics of the massive and the massless 
sectors of the theory (see equations @, (|5|) and @) leads to non-trivial 
jacobian which breaks the BRST invariance [11, 12 1. 

Therefore, due to quantum corrections, the integrability and the BRST- 
invariance properties of the theory fall into two different sectors; the former 
being valid in the massive sector alone and the latter in the meson sector. 
This can be interpreted as the stability of the unphysical (5 particles. We 
have verified this characteristic of the theory in the quasi-classical approx- 
imation by numerical methods. 

The theory simplifies in the large N limit, and meson wavefunctions 
and masses can be computed. Large N corrections are well known in the 
literature, and the numerical problem can be tackled. 

Nevertheless, two dimensional QCD is far from trivial. In spite of the 
large N techniques, the formal aspects of the theory have not been fully 
understood, and only recently the vacuum structure has been studied, and 
separated from the description of the massive excitations. Such structure 
are the core of the understanding of the Schwinger model, and led to very 
profound consequences in that case. Our aim has been to deepen the 
understanding of the theory, obtaining results similar to those known to 
two dimensional QED, the Schwinger model. 

However, it is clear that a development in that direction encounters a 
wall, since the theory is not soluble. Nevertheless, the numerical methods 
used, permit to obtain new structures otherwise unknown. Moreover, there 
are indications that the theory has an unexpected simplification in the 
massless case. 

The result of our paper is to show that in spite of the complexity of 
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the situation, and the fact the the model is used to study QCD mesons, 
the decaying amplitudes are simpler than imagined before, and indicate 
further structures not known before. The zeroes of the decaying ampli- 
tudes are a demonstration of that fact. The theory is not integrable, but 
the decay amplitudes are nearly vanishing. Our numerical results are an 
"experimental" constatation of that fact, and in section () we give a field 
theoretic argument to support that fact. 

Finally, we wish to point out that the methods are not borrowed from 
techniques invented for integrable systems, but rather well defined and 
established techniques based on the computation of the exact fermionic 
determinant, leading to the bosonised version, namely the gauge Wess 
Zumino Witten model and its gauge interaction. 

It is of prime importance to generalise the concept of quasi-integrability 
to higher dimensions. Indeed, Bardeen plj has recently pointed out that 
helicity amplitudes in high-energy QCD are very simple at tree level and 
are described by a self-dual Yang-Mills theory. The classical solution of 
this theory strongly resembles the Bethe Ansatz solution of integrable two- 
dimensional models. Moreover, the one-loop amplitudes are reminiscent of 
those corresponding to anomalous conservation laws. It is known that the 
self-dual Yang- Mills theory is an integrable theory and is described by very 
simple actions p3|. On the other hand, integrable Lagrangians with either 



anomalies [24] or with non-vanishing amplitudes for particle production]25| 
are known and are well documented in the literature. It remains an inter- 
esting open problem to see whether the quasi-integrability idea is the most 
efficient framework for the description of non-trivial dynamics in theories 
with higher conservation laws, in general space-time dimensions, in spite 
of the Coleman-Mandula no-go theorem [26] and its more general version 
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Figure 3: Decay rates of the 5 th , 7 th , 9 th and 11 th excited state into two mesons in 
the 1 st excited state are plotted as functions of the outgoing momenta for massless 
fermions. The vertical bars mark the on-mass-shell values. 
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Figure 4: Decay rates of the 7 th , 8 th , 9 th and 10 th states into two mesons in the 
1 st excited state are plotted for massive fermions. 
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